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Abstract 

Preparation of Schrodinger-cat-like states via conditional output measurement 
on a beam splitter is studied. In the scheme, a mode prepared in a squeezed 
vacuum is mixed with a mode prepared in a Fock state and photocounting is 
performed in one of the output channels of the beam splitter. In this way the 
mode in the other output channel is prepared in a Schrodinger-cat-like state 
that is either a photon-subtracted or a photon-added Jacobi polynomial squeezed 
vacuum state, depending upon the difference between the number of photons in 
the input Fock state and the number of photons in the output Fock state onto 
which it is projected. Two possible photocounting schemes are considered, and 
the problem of monitoring cats that are "hidden" in a statistical mixture of states 
is studied. 



1 Introduction 



The interference of probability amplitudes is one of the most specific features of 
quantum theory and it has been discussed since the early days of quantum me- 
chanics. The famous Schrodinger-cat-like states are a typical example. The cat 
stands for a macroscopic object which may be in a superposition of states corre- 
sponding to macroscopically distinguishable beings (living and dead) PJ. However 
the experimental demonstration of quantum interference effects for macroscopic 
systems is very difficult due to the irreversible interaction of such a system with 
its environment Q. Recent experimental progress has rendered it possible to 
generate Schrodinger-cat-like states on a mesoscopic scale, as it was successfully 
demonstrated in neutron interferometry || and atom optics in a trap Q. 

Despite the large body of work, Schrodinger-cat-like states of optical fields 
have not been observed so far. Recently a sophisticated method for demonstrating 
Schrodinger-cat-like states of traveling optical fields by inferring them from noisy 
data has been proposed ||. The method is based on the experimental scheme 
proposed in || (see also H), and it uses appropriate data processing in order 
to calculate the intrinsic (undistorted) state from the (distorted) state actually 
produced and detected. 

Recently we have shown that Schrodinger-cat-like states can also be pro- 
duced by conditional measurement on a beam splitter. In particular, when a 
mode prepared in a squeezed vacuum is mixed with an ordinary vacuum and a 
(nonzero) photon-number measurement is performed in one of the output chan- 
nels of the beam splitter, then the mode in the other output channel is prepared 
in a Schrodinger-cat-like state ||. The scheme can be extended to the more gen- 
eral case when the mode prepared in the squeezed vacuum is mixed with a mode 
prepared in an arbitrary n-photon Fock state and the photon-number measure- 
ment yields an arbitrary number m. The conditional output states are either 
photon-subtracted (n < m) or photon-added (n > m) Jacobi polynomial squeezed 
vacuum states, i.e., states that are obtained by (\n — m\ times) repeated applica- 
tion of either the photon destruction operator or the photon creation operator, 
respectively, to a Jacobi polynomial squeezed vacuum state ||. It is worth not- 
ing that the two classes of states - similarly to the ordinary photon-subtracted 
(n = 0) U and photon-added (m = 0) ]I0[ squeezed vacuum states - represent 
Schrodinger-cat-like states. 

In this paper we study the properties of these classes of Schrodinger-cat-like 
states, with special emphasis on the experimental conditions. In fact, due to the 
extreme fragility of the quantum interferences the attempt to demonstrate them 
may fail if no rigorous compensation for the detection losses is performed, at least 
in the conditional measurement. We compare photon chopping [[H]] (considered 
in H) with single-detector photocounting (considered in 0) and give an analysis 
of the corresponding data processing algorithms. Further, the effect of losses in 
homodyne detection of the states produced is addressed. 



1 



This paper is organized as follows. In Sec. |2] we outline the basic scheme for 
generating photon-subtracted Jacobi polynomial (PSJP) and photon-added Ja- 
cobi polynomial (PAJP) squeezed vacuum states and briefly address their prop- 
erties. In Sec. H] we compare the two photocounting schemes for conditional 
measurement and analyse the corresponding algorithms used for processing the 
experimental data. Finally, in Sec. f| we give a summary and some concluding 
remarks. 



2 Basis scheme for generating the Schrodinger- 
cat-like states 

The quantum description of the input-output relations of a lossless beam splitter 
are well known to obey the SU(2) Lie algebra | I2| . In the Schrodinger picture, 
the output-state density operator g out can be related to the input-state density 
operator g in as Q ou t = V^ QinV , where V can be given by [TJ] 



y = e -i(<PT-<PR)La e -2i9L 2 e -i{tpT+<PR) L 3 



(1) 



with 
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Let us consider the experimental setup as depicted in Fig. 1. A radiation- field 
mode prepared in a quantum state g m \ is mixed at a beam splitter with another 
mode prepared in a Fock state \n), so that the input-state density operator reads 



as 



Qm{n) = Qua ® |n) 22 (n| 



(3) 



The output-state density operator g out = g out (n) can then be given by 
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From Eq. (f|) we see that the output modes are prepared in a highly entangled 
quantum state in general. When the photon number of the mode in the second 
output channel is measured and m photons are detected, then the mode in the first 
output channel is prepared in a quantum state whose density operator ^ out i (n, m) 
reads 



Qauti(n,m) 
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The probability of such an event is given by 
Tri{ 2 (m|£ out (n)|m) 2 } 

where the abbreviations 

v = n — m, ji = max(0, u), 5 = fi — v (7) 

have been used. 

Let us now consider the case when the first input mode is prepared in a 
squeezed vacuum state £ in i = 5(£)|0)(0|5''''(£), where 



5(O|0> = exp{-i [i{a\f - £*a 2 ] } |0> = (1 - £ ]S ^ ± r «P|2p>, (8) 

£ = |£|e^, k = — e i</3 « tanh |£|. Combining Eqs. (ffl) and (|J) and using Eq. (|[), we 



find that the density operator of the output state reads []10| 

Qouti{n,m) = \ty n>m )(ty n!m \, (9) 

where 

[ alH P(M,«-™)( 2 |T| 2 -l)£(O|0> forn<m, 

*»m) ~ S ~ (10) 

{ (^yP^ h - m \2\T\ 2 -l)S^')\0) foin>m 

is a PSJP (n<m) or a PAJP (n>m) squeezed vacuum state [Pj^(z), Jacobi 
polynomial; f = |£'| e i ( < ^+ 2l ^) ; tanh |«f | = |T| 2 tanh In the photon-number 
basis l^n.m) reads 

Jy "- m £ (k-v)\ \k) 



X 



with = T 2 k. From Eq. (p]) we easily see that when the difference between the 
number n of photons in the second input channel of the beam splitter and the 
number m of photons detected in the second output channel, i.e., the parameter 
v = n — m, is even (odd), then the mode in the first output channel is prepared 
in a PSJP or PAJP squeezed vacuum state \^f n>m ) that contains only Fock states 
with even (odd) numbers of photons. Similarly to ordinary photon-subtracted 
and photon-added squeezed vacuum states [EL H, the PSJP and PAJP squeezed 



vacuum states are Schrodinger-cat-like states. From Eq. ( |TT| ) the normalization 
constant J\f n , m can be calculated to be 

" (-\R\^ fn\" (-\R\y (n\~ (2p + k)\(2p + j)\ 

n ' m h ( k ~»y- w h v - »y w k w&> + »)w |K 1 ' [iZ) 

where e = max(0, [(1 — u)/2]), [x] being the integer part of x. Rewriting Eq. (|l~2|) 

as 

^ m ~\^h \ k )h a-»y- liJw J I h w# 

and using the doubling formula for the gamma function (2p+fc)! = 2 2p+A T(p+^±I) 
r(p + ^ L )/ ^ and the definition of the hypergeometric function, we derive 

■ w\-h <- k -^ wk o-")! Um«'H 1 1 21 2 • 1 w ' 

Note that the hypergeometric function in Eq. flT4] ) is proportional to a Legendre 
polynomial, F(^±±, *±i 1, |/? / | 2 ) = (l-|/? , | 2 )-( fc + 1 )/ 2 P k [(l-|K / |)- 1 /2] ; so that Eq. (|lj) 
may be given by 

» (-|fi| 2 ) fc /n\ A jlhgjj M ^ /j - zA (fe + 2Q!^r; 
x r^LfV^L 3F 2 (Z' + ^, /' + ^, /' + 1, 2Z' - l + 1/2; | K '| 2 ), (15) 



where I' = [(I + l)/2]. The hypergeometric function 3 F 2 (ai, a 2 , 03, &i, 6 2 ; in 
Eq. (|15j) can be calculated using standard routines. In particular, it may be 
expressed in terms of the more familiar hypergeometric function 2 F 1 (a 1 , a 2 , 61; z) 
= F(ai, a 2 , 61; z). The probability P(n, m) of producing the state \^ n , m ) is found 
from Eq. @ to be 

. m!(l-k| 2 ) 1/2 .r , , 

P(n, m) = — \ M n m- 16 

To illustrate the properties of \^? n ,m), m Fig. 2 the quadrature-component 
distributions and the Wigner function of a PSJP squeezed vacuum state for n 
= 1 and m = 4 are plotted [P(n,m) = 3.37%]. From an inspection of the figure, 
the state is seen to exhibit all the typical features of a Schrodinger-cat-like state. 
In particular, a more detailed analysis reveals that the Wigner function is a 
superposition of two quasi- Gaussian lobes with an interference structure between 
them. Since the corresponding formulas are rather lengthy, here we do without 



4 



them and only give an explicit expression of the Husimi function, Q(x,y) 
| (a\ty n>m ) \ 2 /(2n), a = 2~ 1 / 2 (x + iy), which takes the form 
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3 A^-fold photon chopping versus single-detector 
photocounting 

Among the practical problems that may be encountered in an experimental gen- 
eration of the Schrodinger-cat-like states, the losses associated with nonperfect 
photon-number measurement may be of primordial importance. For the sake of 
transparency let us restrict attention to the simplest situation and assume that 
an ordinary vacuum enters the second input of the beam splitter || . In this case 
Eq. (|TTD reduces to 
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[P n (z), Legendre Polynomial]^ and the probability P(0, m) =P(m) of producing 
the state \^ m ) can be given by 
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which is nothing but the probability of detecting m photons in the readout mode. 

At present two types of highly efficient photodetectors are available: the 
linear-response photodiodes suitable for measuring strong signals without single- 
photon resolution and the avalanche photodiodes that may achieve single-event 
discrimination but are then saturated. It has therefore been suggested to spread 
a field whose photon-number statistics is desired to be detected over an array of 
such highly efficient avalanche photodiodes, on using passive optical multiports. 



1 Note that the finite sum in ||] can be expressed in terms of P n (z), on using the relation 
n\J2 l ^ ] (2\x\)- 2k (k!)- 2 {(n - 2fc)!}~ 1 =i n (l - \x\ 2 ) n / 2 \x\- n P n (-i\x\/y/l - \x\ 2 ). 
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Following [[[1], we refer to this method as photon-chopping and to the photodi- 
odes as type I detectors. For a 27V/-port apparatus the probability of recording 
k coincident events when m photons are present is given by 



pM = w [')n-%^-'r pi) 

for k < m, and PN I (k\m) = for k>m. Note that PN I (k\m) — >Sk, m for iVj — >oo. 
In Eq. ([H]) perfect detection is assumed. The effect of nonperfect detection may 
be modelled by placing an absorber in front of the signal before it enters the 2 Al- 
port. This corresponds to a random process such that photons are excluded from 
detection with probability 1—rji, rji being the efficiency of the photodiodes (Note 
that typically rjj = 0.8 . . . 0.94). The probability of recording k coincident events 
then modifies to 

P Nim {k\m) = Y,P*M l ) M iAm), (22) 

I 

where the matrix Mi m (j}i) is given by 

Mi, m {m) = (?W-^r~' (23) 



for l<m, and Mi tJn (j]i) — for l>m. Since detection of k coincident events can 
result from various numbers m of photons, the conditional measurement yields a 
statistical mixture 

Ql{k) = J2PN I>Vl (™\k) \*m){*m\ (24) 
m 

rather than a pure state \^f m ). In Eq. fl24|), P^ IjriI (m\k) is the probability of m 
photons being present under the condition that k coincident events are recorded. 
The conditional probability P^ uVl {m\k) can be obtained using the Bayes rule, 

P NltVl (m\k) = - 1 P Nim (k\m) P(m). (25) 

Here, P(m) is the prior probability fl20D of m photons being present, and accord- 
ingly, P Nl)Vl (k) is the prior probability of recording k coincident events, 

PN im (k) =J2PN 1 , ni (k\m) P(m). (26) 

m 

Let us now consider a so-called type II detector that is able to discriminate 
between zero, one and a few more photons, but with low detection efficiency (rju 
~0.3). Using such a (single) detector for measuring the photon number yields 

Qn{k) =Y. P mi( m \k) |* m )(* m |, (27) 
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where, according to the Bayes rule, the conditional probability P mi {m\k) is now 
given by 

Pvu(m\k) = J J-^M k , m (Vn)P(rn), (28) 

with 

Pm(k) =J2 M k,rn{mi)P{m). (29) 

m 

In order to compare the conditional output states that are produced in the two 
schemes of photon-number measurement, we have calculated the (dimensionless) 
Shannon entropy 

S Kii) = -J2 P i{n)( m \ k ) ^Pi(ii)Hk) (30) 

m 

of the statistical mixtures of states, 

(*) = £ P W"»l*)l*m><ttm|, (31) 
m 

as given by Eqs. ([24]) and (p7j) . The Shannon entropy is a measure of the spread 
of the distribution Pj^^(m\k), i.e., it is a measure of the deviation of Qi^i) from 
a pure state. Note that for a pure state SV(77) = is valid. From Fig. 3 we see 
that 5*7- can always be reduced below Su when the number of type I detectors 
in the photon chopping scheme is sufficiently increased. Hence photon chopping 
(with type I detectors) may be much more suitable for preserving the quantum 
interference features in the conditional output state than the use of a single type 
77 detector. Photon chopping yields a mixed state that is less spread and "more 
pure" than the state obtained by using a single type II detector for photon- 
number measurement. Clearly, when in the photon chopping scheme type // 
detectors are used in order to record coincident events, then this scheme is less 
suitable for photon counting than a single type // detector. For r/j = r\ n the 
two schemes yield equal conditional output states only in the limit when Nj — > 
oo. For finite Nj there is always a nonvanishing probability that the number of 
recorded coincident events is smaller than the number of photons, Eq. (pT|) . In 
particular with regard to pure-state generation we have 



lim lim Si 



0. 



lim Sji 

mi-* 1 



0. 



(32) 



The conditional output states (|24]) and (|27D can be determined using balanced 
homodyne detection and measuring the quadrature-component distributions 



Pi(ii)(x,<p\k) =J2Pi(n)(m\k)p(x,(f\m), 



(33) 
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where |J 

p(x,(p\m) = \(x,ip\^ m )\ 2 
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with A = 1 + |k / | 2 + 2|k / | cos(2(/) — </v) and X - 
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)/A . In Fig. 4 



we report the results of simulated measurements. In the case of a single type 
// detector, Fig. 4(a), the quantum interferences are totally smeared and non- 
observable. Although somewhat smeared, the quantum interferences are observ- 
able in a photon-chopping scheme with type / detectors, Fig. 4(b). The result 
obviously reflects the above mentioned fact that the conditional output state Qu, 
Eq. fl2"T|), is "more mixed" than the state qi, Eq. (|2"4]), in general. 

The reconstruction of the states |\l/ m ) from the mixed state @n{k), Eq. 
can be achieved using the inverse Bernoulli transform fl5|, 




Similarly, inverting Eq. (|24]) , we obtain (iVj — >oo 



(35) 
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where (P, 



Ni)m,k 



is the inverse of the matrix (P, 
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lowing recursion relation is valid [111 
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(37) 



Note that Eq. (j36|) reproduces the exact components |\l/ m )(\I/ m | only for Nj ^oo. 

Applying Eq. (|33| ) and Eq. (^) (for finite Nj), the states |\l/ m ) can be recon- 
structed from the homodyne data of the measured mixed states. Examples of 
reconstructed quadrature-component distributions are shown in Fig. 5. From the 
figure we see that processing the homodyne data according to Eqs. ( p5|) and (|36|) 
allows us to restore the quantum interferences in the quadrature-component dis- 
tributions of the component states of the produced statistical mixtures of states. 
It should be noted that with increasing number of measurements the Bernoulli 
inversion in Eq. ( |35"D yields the almost perfect interference structure [Fig. 5(a)]. In 
order to realize (for the same number of measurements) a comparable accuracy on 
the basis of Eq. (|36|), a sufficiently large number Nj of channels (type /detectors) 
in the photon-chopping scheme must be used [Nj = 50 in place of Nj = 20 used in 
Fig. 5(b)]. This is obviously due to the fact that in the photon-chopping scheme 
the probability that a photon impinges on an already saturated photodiode ap- 
proaches zero only for Nj — ► oo. On the other hand, photon chopping already 
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yields reasonable results for small amounts of data, even when the number of 
channels is reduced [Fig. 5(d)]. From Fig. 5(c) it is seen that the error in the 
single- (type-//-) detector scheme drastically increases with decreasing number of 
measurements. This result tells us that photon-chopping may be more powerful 
than a single-detector scheme when the produced state tends to a macroscopical 
one and the amount of data needed becomes great. 

Finally, let us mention that with respect to the probability of producing par- 
ticularly macroscopic Schrodinger-cat-like states the chopping scheme may be 
more suitable than the single-detector scheme. In fact, from Fig. 6 we see that 
the probability of recording k > clicks is always higher for the photon-chopping 
method as for the single-detector scheme. This effect can be made a bit more 
pronounced choosing larger Nj. 

Let us briefly comment on the use of a nonperfect homodyne detector for 
measuring the quadrature-component distributions (^). Since in a realistic ho- 
modyne experiment the quadrature components cannot be measured exactly, we 
may assume that instead of p/(//)(x, <p\h) smeared distributions 

p I{II) (x,(p;r]\k) =J2 p i(ii)( m \ k )p( x i ( P',v\ m ) (38) 

m 

are measured, where 

/oo 
dy f(x-y,T])p(y,ip\m), (39) 
-oo 

f(x; rj) being some positive single-peaked function of x, such as a Gaussian, 



(r/, quantum efficiency of the homodyne detector, with 77 < 1). Combining Eqs. (|3 
), and ( |4"0| ) and performing the ^/-integration yields 
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(41) 



p = (1 - \k'\ 2 )/A + 1/(2(7). It can be readily seen from Eqs. @ - (|T]) that 
p(x, <p; r]\m) — > p(x,(p\m) for 77 — > 1. The notorious fragility of the interference 
structure may be directly seen by comparing the quadrature-component distri- 



butions fl34|) and (41). For example, for m = 3 (an a mean number of photons (n) 
= 15) the interference structure is completely smeared out if r] < 0.94 (Fig. 7). 

It is well known that for compensating the losses of the homodyne detector 
one can make a detour via the density matrix in the Fock representation. The 
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density-matrix elements are reconstructed from the noisy quadrature- component 
distributions using loss- compensating kernels, 



Qn,n'{m) = (n\^ m ){9 m \n') = / dip dxp(x, ip\ rj\m)K njn >(x, ip; rj), (42) 



K n>n /(x, ip; rj) being given in |L3] (note that such a compensation is possible only 
for 77 > 0.5). Alternatively one can first reconstruct the density-matrix elements 
in the Fock basis using the kernel functions for perfect detection and then apply 
an inverse Bernoulli transform to reconstruct the true density-matrix elements 



We have shown that quantum-state preparation via conditional output measure- 
ment on a beam splitter can be advantageously used for preparation of a great 
variety of Schrodinger-cat-like states. When a mode prepared in a squeezed vac- 
uum and a mode prepared in an arbitrary Fock state are superimposed by a beam 
splitter and an arbitrarily chosen number of photons is recorded in one of the 
output channels of the beam splitter, then the mode in the other output channel 
is prepared in either a photon-subtracted or a photon-added Jacobi polynomial 
squeezed vacuum state, the latter being obtained by applying an operator-valued 
Jacobi polynomial to a squeezed vacuum state. All the PSJP and PAJP states 
represent examples of Schrodinger-cat-like states, provided that the sum of the 
number of incident photons and the number of detected photons is nonzero. 

In order to project onto Fock states, we have studied two methods of di- 
rect photon counting which may be realized with currently available techniques: 
single-detector photon counting and iV-fold photon chopping. Both methods pro- 
duce statistical mixtures of Schrodinger-cat-like states rather than pure states in 
general, because of nonperfect detection. We found that photon chopping of- 
fers the possibility of direct observation of the quantum interferences. Moreover 
this method can be advantageously used to reduce the amount of data needed 
for reconstructing the Schrodinger-cat-like states in the produced mixed state, 
which can be measured by balanced homodyning. When the number of recorded 
data is suitably large, then the use of a single (low-efficiency) detector may be 
more advantageous. In this case the pure-state components of the statistical 
mixture of states can simply be calculated from the measured data using the 
inverse Bernoulli transform. To realize the same accuracy in photon chopping, 
the number of channels and detectors must be relatively high. 
Acknowledgements This work was supported by the Deutsche Forschungsge- 
meinschaft . 
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Figure 1: Experimental setup for generating the Schrodinger-cat-like states. A 
signal mode prepared in a state Qmi = S(£) | 0) (0|S^(£) is mixed at the beam splitter 
BS with another input mode prepared in a Fock state |n), and m photons are 
recorded by the detector D in one of the output channels of the beam splitter. 
The quantum state g out i of the mode in the other output channel is found to be 
"collapsed" to either a PS JP {n <m) or a PAJP {n > m) squeezed vacuum state. 
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Figure 2: Quadrature distribution (a) and the Wigner function (b) of a PSJP 
squeezed vacuum \^ n>m ) [n—1 and m = 4] for k' — -0.81 [|T| 2 = 0.9, |/t|=0.9]. 
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Figure 3: The Shannon entropies Sj and Sjj of the mixed states (|24]) and (|27j). 
respectively. In (a) it is assumed that k = 3 events are recorded and the quantum 
efficiencies are 77/ = 0.85 and 77/7- = 0.3. In (b) the parameters are k = 5 and 77/ = 0.5 
and r\n = 0.3. The calculations are performed for k' — — 0.7 (|T| 2 = 0.9, \k\ =0.77). 
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Figure 4: The quadrature-component distributions pn(x, <p\k) (a) and pi(x, ip\k) 
(b), respectively, of the mixed states © (77/7 = 0.3) and (||) (iV 7 = 20, 771 = 0.9) 
are shown for ip — and /c = 3. The theoretical distributions (|34] ) (solid lines) 
are compared with the histograms obtained from computer-simulations of 10 6 
measurements. The calculations are performed for /t' = —0.81 (|T| 2 = 0.9, \k\ = 
0.9). 
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Figure 5: The quadrature-component distributions p(x,(p\m) of \^Sf m ) recon- 
structed from the homodyne data, (a,c) according to Eq. (|35| ) (rju = 0.3) and 
(b,d) according to Eq. (|36|) (Nj = 20, T]j = 0.9), are shown for <p = Q and m = 3. 
The theoretical distribution fl3HD (solid line) is compared with the distributions 
reconstructed from (a,b) 5 • 10 5 measurement, (c) 10 4 measurements, and (d) 10 3 
measurements. The calculations are performed for k' — —0.81 (|T| 2 = 0.9, |«| = 
0.9). 



16 



0.8 



CO 

•2 0.6 



42 
O 

ft 0.4 

_o 



0.2 



2 
A: 



Figure 6: The prior probabilities fl26|) (solid bars) and fl29[) (dashed bars) for 
recording k clicks in the photon-chopping scheme = 0.9, Nj = 20) and the 
single-detector scheme (r/jj = 0.3), respectively, are shown. The calculations are 
performed for k' = -0.81 (|T| 2 = 0.9, \k\ = 0.9). Note that for k = 3 the type i7 
detector only gives 12% of the probability obtained by photon chopping. 
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Figure 7: The theoretical quadrature-component distribution (|34j) at the phase 
<p = (solid line) is compared with the quadrature-component distribution ([41]) 
obtained in nonperfect homodyne detection for 77 = 0.98 (dashed line) and rj = 0.94 
(dotted line). The values of the other parameters are as in Fig. 3 The calculations 
are performed for k' = —0.81 (|T| 2 = 0.9, |k|=0.9). 
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